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It is shown that the electronic conduction in silicon-on-insulator (SOI) layers exhibits a metallic 
regime which is very similar to that in high- mobility Si- metal oxide semiconductor structures (MOS). 
The peak in the electron mobility versus density, the strong drop in resistivity and the critical 
concentration for the metal-insulator transition are all consistent. On the basis of our SOI data for 
the temperature and in-plane magnetic field dependence of the resistivity, we discuss several models 
for the metallic state in two dimensions. We find that the observed behavior can be well described 
by the theory on the interaction corrections in the ballistic regime. For the investigated regime, the 
temperature dependent screening of scattering potentials gives also a good description of the data. 



I. INTRODUCTION 



The discovery of an apparent metallic state in Si- 
metal oxide semiconductor (MOS) structures!!] has at- 
tracted much attention as it seemed to contradict the 
two-dimensional (2D) localization behaviom. It was 
found that above a critical carrier concentration n c , the 
resistivity strongly decreases with decreasing tempera- 
ture T (metallic regime), whereas it increases for lower 
densities n < n c (insulating regime). The scaling param- 
eter To shows a critical behavior around njp. 

At first the discovery of the metallic state in two di- 
mensions was met with great scepticism. But the metal- 
to-insulator transition (MIT) was confirmed first in Si- 
MOS structures from a different sourceo and further in 
several other 2D material systems, which are n- and p- 
SiGe and AlGaAs and n-AlAs structures (for references 
seeo). 

After these experimental findings, the question was 
raised whether the metallic behavior manifests a new 
electronic ground state with dominating quantum behav- 
ior or if the strong resistivity drop towards low T is based 
on familiar_(i.e. non coherent) effects (for a detailed dis- 
cussion seeo and references therein). 

In this work, we report on the finding of a metallic state 
in silicon-on-insulator (SOI) structures. In section |l[ on 
the experiment, we describe the investigated samples and 
the setup. We further present the data on the tempera- 
ture and the magnetic field dependence of the resistivity 
and we show the mobility versus density behavior. The 
experimental results of the SOI system are compared in 
detail with the properties of Si-MOS structures and the 
differences are discussed. The discussion of the results 



follows in section III. Here we discuss and compare five 



important suggestions on the physical origin of the metal- 
lic state in two dimensions. It turns out that the co- 
herent interaction corrections to the conductivity in the 
ballistic regimeBB are able to describe the observed re- 
sistivity behavior quite well. For the investigated regime 
also the temperature dependent screening of scattering 
potentials!!! gives a fairly good description of the data. 
After a short summary, we describe in appendix [X] our 
numerical calculations on the charged trap state model 
in detail. 



II. EXPERIMENT 

A. Samples and setup 

Our investigations were performed on two high- 
mobility SOI structures which were prepared recently in 
Finland. The SOI Hall bars have identical peak mobil- 
ities of /ip = 11, 600 cm 2 /Vs. They were fabricated on 
commercial 'Unibond' wafers produced by the Smart-Cut 
process!! The thickness of the Si-layer amounts 200 nm 
and that of the buried oxide layer 400 nm in these wafers. 
The background doping is n-type of about 10 15 cm~ 3 
which freezes out at low temperature. After a two-step 
thermal oxidation, the final thickness of the Si-layer is 
60 nm. The gate oxide thickness is about 90 nm with a 
dielectric strength of approximately 1 V/nm. The electric 
contacts were phosphorus-implanted to an activated car- 
rier concentration of 4 x 10 19 cm~ 3 . The Al-metallization 
is 300 nm thick and sintered at 400°C. A scheme of the 
final layer sequence is shown in Figure |. The etched Hall 
bars have a size of 800 /im times 100 (im giving a length 
to width ratio of 8. All resistivity and Hall measure- 



corresponding author: Tel.: +43-732-24689603; Fax.: +43-732-2468-8650; e-mail: G.Brunthaler@hlphys.uni-linz.ac.at 



1 



ments were performed in a four terminal AC-technique 
at low frequencies of typically 13 or 17 Hz. The SOI struc- 
tures were investigated in a 4 He cryostat down to 1.4 K. 
Most investigations were performed on sample number 
L2, sample LI showed in several measurements practi- 
cally the same behavior. 

The SOI sample number L2 was intensively inves- 
tigated in a density range from 1.3 x 10 11 to 1.04 x 
10 12 cm -2 . At the low density side of this range, the 
contacts became high ohmic. Due to capacitive effects, 
the measured voltages showed a delay time of more than 
100 ms before they came into a steady state. Thus voltage 
measurements in DC-technique lead to other values than 
those with the 13 Hz AC-technique. We limited our mea- 
surements to the range were the AC- and DC-tcchniquc 
gave the same results and performed all investigations in 
the AC-technique as this gives a better signal-to-noise ra- 
tio. At the higher limit of the investigated density range, 
the gate started to leak and the current between gate and 
sample would have lead to errors in the resistivity calcu- 
lations. The investigations were thus kept in the above 
given density range. 
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FIG. 1. Scheme of the sample structure for sili- 
con- on- insulat or structures. The vertical axis shows the layer 
sequence (not to scale), the horizontal axis the band offsets. 



B. Resistivity behavior 

From the resistivity and Hall measurements, the car- 
rier density n and the mobility fi were calculated assum- 
ing the linear Drude behavior to be valid. This is not 
clear a priori as quantum interference effects might give 
relevant contributions to the conductivity and would al- 
ter the necessary evaluation. The possible importance of 
quantum effects will be discussed later in this work. For 
now, n and /i should be used as apparent Drude param- 
eters for the 2D electron gas. 

We have evaluated the dependence of mobility /i on 
the electron density n by varying the gate voltage V g 
at a constant temperature of T = 1.5 K. In Fig. 0, the 
/i(n)-dependence forJSOI is compared with that for sev- 
eral Si-MOS sampled^ with different mobilities. The SOI 



structure has a peak mobility \x v of 1 1,600 cm 2 /Vs at a 
density of 6.5 x 10 11 cm~ 2 . As can be seen from Fig. ||, 
the Si-MOS samples have lower or higher peak mobil- 
ities, depending on the sample quality, but the overall 
behavior of fJ,(n) is very similar for the two sample types. 

At high densities, the electron wave function is 
squeezed by the strong electric field of the triangular 
potential towards the gate-sided Si/Si02 interface and 
thus interface roughness scattering dominates. At low n, 
the dominating scattering process is caused by impuri- 
ties. At the transition region between the two scattering 
mechanisms, the mobility reaches its peak value /Xpli9. 




FIG. 2. Mobility fj, versus density n for a SOI structure in 
comparison with several Si-MOS samples at a temperature 
of 1.5 K. 

The most typical feature of the metallic state in 2D 
systems is expressed by the strong resistivity drop to- 
wards lower temperatures. This strong drop in p(T) is 
clearly visible in Fig. ||b for the SOI sample in the range 
of 1.67 x 10 11 < n < 1.04 x 10 12 cm~ 2 . The maximum 
decrease in p(T) amounts up to about a factor 3.5 for 
n = 2.47 x 10 11 cm~ 2 . For the two lower density curves 
with n = 1.28 and 1.44 x 10 11 cm -2 , an insulating behav- 
ior is observed. The critical concentration n c , at which 
the behavior changes from insulating to metallic, lies at 
about 1.55 x 10 n c m - 2 . 

For comparison, in Fig. ||a and the very high and 
medium mobility Si-MOS samplesEll Si-15 and Si-4/32 
with Hp = 27,000 and 8, 000 cm 2 /Vs, respectively, are 
shown. In sample Si-15, the strong decrease in p(T) is 
shifted to lower T and the maximum decrease amounts 
up to a factor 7, in Si-4/32 the decrease is at somewhat 
higher T in comparison to the SOI structure and the 
maximum decrease is about a factor 3. This shows a 
clear trend of the properties of the metallic regime with 
the sample quality, i.e. the peak mobility fip an d mani- 
fests the similarities between SOI and Si-MOS structures. 
Additionally, in Fig. dashed-dotted lines at Ej — ksT, 
AksT and 16fcsT are shown. 
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FIG. 3. Resistivity versus temperature behavior in units of h/e 2 for a) high-mobility Si-MOS sample Si-15 
Op = 2.7m 2 /Vs), b) SOI-L2 sample (1.17m 2 /Vs), and c) medium mobility Si-MOS Si-4/32 (0.8m 2 /Vs). The 
electron densities for Si-15 are between 0.93 and 10.2 x 10 11 cm" 2 , for SOI-L2 are 1.28, 1.44, 1.67, 2.02, 2.47, 3.26, 
4.06, 5.65, 7.50 and 10.4 x 10 11 cm" 2 and for Si-4/32 between 2.34 and 28.9 x 10 11 cm -2 . 



Figure |J shows a direct comparison of the strong 
changes in p(T) for the four samples Si-15, Si-43, SOF 
L2 and Si-4/32 for the same low-T resistivity p — 
0.0274/i/e 2 . The same behavior, as described above, is 
visible: the higher the peak mobility, the stronger is the 
change in p(T) and the lower is the temperature where 
the increase starts. 
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FIG. 4. p(T) for samples Si-15, Si-43, SOI-L2 and Si-4/32 
for one and the same low-T resistivity of 0.0274/i/e 2 . The 
crosses on the curve from left to right mark the temperature 
where k B T = Bp/16, E F /4, and E F . 

For Si-MOS structures, it was shown that the critical 
conductance g c = l/p c changes systematically with the 



peak momentum relaxation time fp, = m*p p /e, which 
is a measure of the sample qualityli3. Here m* denotes 
the effective conductivity mass which is the transverse 
mass of the two lower lying valleys of ml — 0.19 m e and 
e the elementary charge. For the SOI structure, g c is 
about l.le 2 /h and is compared_m Fig. || with the val- 
ues for several Si-MOS samplestx It can be seen, that 
the g c -value for SOI coincides quite well with the general 
behavior of the Si-MOS samples. 
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FIG. 5. Critical conductance g c in units of e 2 /h versus 
the inverse peak momentum relaxation time l/r p for the 
SOI structure L2 (star symbol) and several Si-MOS samples 
(square symbols). The dashed line shows a linear regression 
curve. 
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As was seen so far, the SOI structure behaves very 
similar to a high mobility Si-MOS sample. Therefore, we 
ask the question: What are the similarities and differ- 
ences between the two sample types? Due to the different 
effective masses to* perpendicular to the Si/Si02 inter- 
faces and due to strain effects, the six-fold conduction 
band (CB) degeneracy at the A-points of bulk Si is lifted 
and only the two lower lying valleys with their longitudi- 
nal axis perpendicular to the 2D layer are occupied. The 
Si layer in SOI is only 60 nm thick, below there is the 
buried oxide layer, which constitutes also of a thermally 
grown oxide fabricated before the wafer bonding process. 
The active Si-layer has an n-type background doping of 
about I0 15 cm~ 3 which freezes out at low T. In the Si- 
MOS samples, the electrons are confined in the surface 
inversion layer of a nominally p-type doped bulk Si crys- 
tal. Our high mobility Si-MOS samples were doped with 
an acceptor density of about 2 x 10 15 cm -3 . 

Due to the applied gate voltage V g , an approximately 
triangular potential well results for the electrons at the 
upper Si/Si02 interface. Nevertheless, the extension of 
the electronic wave function is different for the SOI and 
the Si-MOS structure due to the different background 
doping. In the n-type SOI structure, the Fermi level Ef 
is pinned underneath the 2D layer between the donor 
level and the CB edge. The band bending below the 2D 
layer is therefore quite small. In p-type Si-MOS, there 
is a wide depletion layer of about 0.85 /im below the 2D 
layer, as Ep is pinned here between the acceptor level 
and the valence band (VB) edge in the undistorted p- 
type region. The different band bending below the 2D 
layer has an influence on the extension of the wave func- 
tions especially for small electron densities. 
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FIG. 6. Width of electron wavefunction versus electron 
density n for SOI (circle symbol) and Si-MOS (triangle). 
The width is defined to include the center 90 percent expec- 
tation probability of the wave function, 5 percent at the left 
and the right hand side are omitted. 

We have calculated the wave function ^>(z) and the 
triangular potential V(z) by solving the Schrodinger and 



Poisson equations in a self consistent way (with z the co- 
ordinate perpendicular to the interface layer). In order 
to solve the Schrodinger equation, the charge distribu- 
tion p(z) (x \^f 2 (z)\ was taken into account. No exchange 
or correlation effects among the electrons were consid- 
ered, thus giving a first order approximation of ^(z) and 
V(z). Figure ^ shows the width of the wave function ver- 
sus electron density n. The width of \1/ is defined as the 
length where the central 90 percent of the expectation 
probability \^f 2 \ is contained, with 5 percent at the left 
and at the right side ignored. For a high electron density 
of 10 12 cm -2 the width is with 7.6 and 9.2 nm nearly the 
same for both structures. 

At lower densities, the difference in the width of 
increases. The depletion width of 0.85 /im at an accep- 
tor density of 2 x I0 15 cm~ 3 corresponds to a 2D charge 
density of 1.7 x I0 n cm~ 2 . Thus, at that 2D electron 
density, the electric field at the Si/Si02 interface is in 
Si-MOS twice as high as in the SOI structure. As can 
be seen from Fig. ^|, the difference in the wave function 
extension is indeed large at that density. At still lower 
density, the width saturates in Si-MOS, as the field from 
the conducting 2D electrons is less important than that 
from the acceptors in the depletion layer. On the oppo- 
site, the width in the SOI structure further increases, as 
there is no electric field confinement from other charges 
than the conducting electrons itself. 

Figure ^ shows for comparison the potential shape 
V(z) and the wave function ^f(z) at n = 10 11 cm~ 2 . At 
that density, the confining potential is much steeper in 
Si-MOS and thus the energy Eq of the lowest electronic 
subband is with 31 meV much higher than in SOI with 
14 meV. In Fig. |7|, the difference in the width of the wave 
functions is also clearly visible and is caused by the differ- 
ence in the potential steepness. For the case that higher 
subbands ^ n (z) are occupied, the differences in their en- 
ergies and extensions should be even larger. 
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FIG. 7. Comparison of the wave functions ^(z) and the 
potential shape V(z) for the SOI (solid line) and Si-MOS 
structure (dashed line) at an electron density of lx f n cm~ 2 . 
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In the considered density range, the width of the wave 
function is always smaller than the Si-layer thickness 
of the SOI structure and thus the electron do not feel 
the lower Si/SiC>2 interface towards the buried oxide 
layer. In order that such an effect could become signifi- 
cant, the Si layer would have to be as thin as 20 nm for 
n w 1 X 10 11 cm~ 2 and even thinner for higher densities. 

C. Magnetoresistivity 

We have also investigated the SOI structures in per- 
pendicular and parallel (i.e. in-plane) magnetic field B. 
Figure [sj shows typical magnetoresistivity curves p xx (Bj_) 
for B perpendicular to the 2D layer for different densities 
at T — 1.4 K. For the lower densities, the weak localiza- 
tion peak around B — is clearly visible. At about 
1 T, the Shubnikov - de Haas oscillations set in and gain 
strongly in amplitude towards higher magnetic fields. At 
around 4 T, the spin splitting is visible. The behavior in 
perpendicular magnetic field is again very similar to that 
of Si-MOS samples with comparable mobility. 



electrons^ and later io the increase in disorder caused 
by the magnetic fielda A strong change of p{Bu) was 
also observed in p-type GaAs/AlGaAs structures £3 
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FIG. 9. Relative change of the magnetoresistivity 
p xx (B^)/ p xx (0) for magnetic fields oriented parallel to the 
2D layer for n = 1.28, 1.44, 1.67, 2.02, 2.47, 3.26, 4.06, 5.65, 
7.50, and 10.4 x 10 11 cm -2 (from top to bottom) at a tem- 
perature of 1.4 K. 
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FIG. 8. Magnetoresistivity p xx for magnetic fields ori- 
ented perpendicular to the 2D layer for n = 3.26, 4.06, 5.65, 
7.50, and 10.4 x 10 11 cm -2 at a temperature of 1.4 K. 

For in-plane magnetic field, the relative change 
Pxx(B\\)/p xx (0) of the magnetoresistivity of SOI sample 
L2 is shown in Fig. || for different densities at T = 1.4 K. 
At high densities of about 1 x 10 12 cm~ 2 , the changes 
are quite small, whereas for the lower densities there is a 
strong increase in p xx by up to a factor 4. At the lowest 
density of 1.28 x 10 11 cm -2 , the observed change is again 
lower than for the previous densities as can be seen from 
the crossing of the upper curves in Fig. |^. This indicates 
a saturation of the strong increase in p xx (B») at about 
1.5 x 10 u cm- 2 . 

An even stronger p(_B||)-dependence was observed 
in higiujpapbility Si-MOS structures in the metallic 
regimeEjo. Again the changes were stronger for the 
lower densities n. The increase with Bu was first 
attributed to the spin polarization of the conducting 



III. DISCUSSION 

Above, we have presented new results on the 2D metal- 
lic state in silicon-on-insulator structures and compared 
its properties with that of high- mobility Si-MOS samples. 
As the physical origin of the 2D metallic state is still an 
open question, we like to discuss the different possibilities 
in detail. 

There appear to remain five important suggestions for 
the 2D metallic state at low temperatures. These are (i) 
a new quantum mechanical ground state for electrons in 
a solid, (ii) scattering of electrons at charged hole traps 
in the oxide layer, (iii) interband scattering between dif- 
ferent bands, (iv) temperature dependent screening of 
impurity scattering, and (v) coherent interaction correc- 
tions in the ballistic regime. If the first suggestion proves 
to be correct, the 2D MIT would indeed be the manifes- 
tation of some new physics, whereas the next three sug- 
gestions are based on settled physical phenomena, which 
may be described by terms like conventional, semiclas- 
sical or non-coherent (on the first sight). The fifth sug- 
gestion was proposed quite recently and it has to be dis- 
cussed whether it can be classified in the terminology 
used above or not. 



A. New quantum mechanical ground state 

The first suggestion (i) of a new quantum mechani- 
cal ground state is based on the scalability of the re- 
sistance in the neighborhood of the apparent MIT in 
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high-mobility Si-MOS structures!]. It was found that 
the p(T, n)-curves for n < n c and for n > n c can be 
brought to coincidence on two distinct curves for the in- 
sulating and the metallic regime, respectively, by scaling 
the temperature with T/Tq(h). The scaling parameter 
exhibits a power law To oc \S n \ b , with 8 n = (n — n c )/n c 
and b ~ 1.6. Such a critical scaling behavior indicates 
the metal- insulator transition to be a quantum phase 
transitior£3o. 

The seminal paper on the absence of quantum diffu- 
sion in two dimensions! 2 ] was based on the assumption of 
non-interacting electrons. In this case the weak local- 
ization contribution decreases the conductivity and leads 
to continuing higher resistivities with decreasing temper- 
atures. In the recent quantum-phase-transition models 
for the unexpected MIT in 2D, the jfilectron-electron (e- 
e) interaction plays a dominant rold_J~EHI. The quantum 
coherent interaction corrections consist of positive and 
negative contributions to the conductivity and may lead 
to a decrease of the resistivity under certain conditionsEil. 
In any case, quantum coherence is a prerequisite that the 
MIT can be a quantum phase transition. 

A metallic behavior, i.e. a strong drop in the resis- 
tivity, was not only observed in the near vicinity of the 
critical concentration n c but also far away, for densities 
up to rt > 1 x 10 12 cm~ 2 (i.e. more than 10 times n c ). As 
the decrease in p(T) looks quite similar near n c and far 
away, it is quite instructive to look for quantum .effects 
over an extended density range. In a recent workEil, the 
temperature dependence of the weak localization was in- 
vestigated in the metallic state of high-mobility Si-MOS 
structures in order to find the borders of quantum co- 
herence. In that work, the temperature threshold T q for 
the single-electron quantum coherence was deduced from 
the crossing point of the temperature dependence of the 
phase coherence time t v and the momentum relaxation 
time r. When t v > r, the coherence time is long enough 
that an electron can return in the coherent state to its 
origin whereas for t v < t this is not noasible. 

It was found for the Si-MOS systemEJ that for n > 2 x 
10 11 cm -2 , the strong decrease in p(T) takes place above 
T q and thus should not be related to a single-electron co- 
herent effect. Also the second border /esT 0C — H/t lies 
below the strong resistivity drop at still higher densities. 
Therefore also the e-e interaction induced quantum cor- 
rections to conductivity should not be responsible for the 
resistivity drop in this regime. This should at least be de- 
duced in the picture where quantum coherent effects and 
screening are attributed to two different physical phe- 
nomena. But^ecently an interesting paper by Zala et al. 
has appearedB, which may indicate that the two before 
mentioned effects are caused by the same physics. This 
will be discussed in more detail later in this work. 

In any .-case, the work on the borders for quantum 
coherence^ is strictly applicable only in the higher den- 
sity regime and does not give a direct answer on the im- 
portance of quantum coherence in the close vicinity of 
the critical density n c . In a recent work, Punnoose and 



Finkelstein showed by a renormalization group analysis 
that especially for high-mobility Si-MOS structures, the 
presence of two degenerate conduction band valleys, sta- 
bilizes and enhances coherent e-e interaction effects in the 
vicinity of n c . They are able to describe the temperature 
dependent changes in conductance in an n-range between 
0.83 and 0.94 x 10 11 cm -2 without any adjustable param- 
eter. So it is quite surprising that at low n the whole drop 
in p seems to come from quantum coherent interaction 
effects while at high n the drop takes place without quan- 
tum coherence and nevertheless the strong T-dependcncc 
looks so similar in the different regions. It is not yet clear 
whether the metallic and insulating behavior belong to a 
quantum phase transition or not. This work focuses on 
the metallic regime and will not be able to judge on the 
metal-to-insulator transition itself. 



B. Charged hole traps 

The second suggestion (ii) on the scattering of elec- 
trons at charged hole trraps in the oxide layer by Alt- 
shuler and Maslov (AM)E3 is prepared especially for Si- 
MOS structures with the Si/Si02 interface on the gated 
side of the 2D layer. The same model should be applica- 
ble also for our SOI structures, as an equivalent Si/Si02 
interface is present for the 2D electron confinement. 

It is well known that in the SiOaJnterface layer there 
exist several types of defects statesEj. In the AM model, 
it is assumed that a relative large number of hole trap 
states exists at a certain trap energy E t , which is spa- 
tially homogeneously distributed in the oxide layer. It 
is further assumed that these trap state lie at a well de- 
fined energy E t o if no external electric field is applied. 
The potential gradient due to an applied gate voltage V g 
causes a linear increase of the trap energy position. In 
order to describe the total energy of the trap states, also 
the mirror charge of the traps inside the 2D electron gas 
has to be included. This leads to a down bending or 
the energetic position towards the interface and causes 
a maximum in the total trap energy E t {z), with z the 
position inside the oxide layer (see Fig. [l8| in Appendix 
[A|). When the Fermi energy lies below certain trap states 
at energies E t (z) (for low n) these states are positively 
charged and cause a large scattering probability for the 
2D electrons. At high densities all trap states lie below 
Ep and only due to the Boltzmann occupation tail, an 
exponentially small number of traps contributes to the 
electron scattering. The charged trap states form dipoles 
together with their mirror charge in the 2D layer which 
has to be taken into account in the calculation of the 
scattering cross section. 

In order to calculate the resistivity, one needs to know 
how the Fermi energy Ep varies with the temperature T 
(in the work of AM, the chemical potential at T = is 
called the Fermi energy Ep whereas the chemical poten- 
tial at T > is denoted by p{T)\ in this work we use 
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Ep(T) instead). AM suggest two scenarios for the T de- 
pendence of the Fermi energy of the 2D electron gas : (a) 
it coincides with the Si substrate or (b) it is decoupled 
from the substrate. 

For case (a), the calculation of p{T) by AM show a be- 
havior which looks quite similar to the observed tempera- 
ture dependence near the critical density n c for high mo- 
bility Si-MOS structures. The strong increase in p with 
increasing T occurs around fc^T w Ep — E t as then the 
occupied trap states are partially emptied (and charged) 
due to the broadening of the Fermi occupation function 
and the scattering rate is increased accordingly. For case 
(a), the insulating and metallic regions are distributed 
quite symmetrically around the separatrix which occurs 
for Ef = max(Et(z)). 

For case (b), the resistivity is decreasing in both, 
regimes (Ep ^ E t ) according to the calculations of AME2I 
and seems to be unable to explain the experimentally 
observed behavior. As mentioned by AM, the more re- 
alistic case is (b) as the depletion layer isolates the 2D 
electron gas from the substrate and both regions should 
have different ouasi Fermi energies. This is also verified 
experimentally!^!. This situation causes some difficulty in 
order to apply the charge trap model to the high-mobility 
Si-MOS structures as well as to the SOI samples consid- 
ered in this work. 

We have performed numerical calculations of the AM- 
trap model for the SOI data. These calculations are de- 
scribed in appendix A in detail. If one applies an in- 
creasing gate voltage V g in order to increase the elec- 
tron density n in the 2D layer, the spatial maximum in 
the trap energy shifts quite fast to lower energies. This 
leads to a strong occupation of the trap states as soon 
as the maximum shifts below the Fermi energy Ep and 
the scattering efficiency decreases drastically. In Fig. 
we have calculated the p(T) dependence at an interme- 
diate temperature of Tj = 3.1 K, where we have a set 
of measurement data. The temperature 7$ was chosen 
so that p(n,Ti) is for a large density interval neither in 
the low-T saturation nor in the "normal" high-T regime. 
The resistivity behavior should be dominated by thermal 
excitation of carriers according to the AM-trap model. 
The low-T saturation of p(T) in the metallic regime can- 
not be described with the AM-trap model alone, as all 
traps get filled according to the exponential Boltzmann- 
tail and the resistivity tends to zero. Only by assuming 
additional residual scattering centers, the resistivity sat- 
urates. We have evaluated the measurement data at Tj 
only for densities where the low-T saturation should not 
have a large influence yet. 

Figure [f(] shows that the decrease in resistivity vs. 
density is too strong by many orders of magnitude in 
the AM-trap model. At n = 2.5 x 10 n cm~ 2 , p should 
be below 10~ 30 h/e 2 . This low value results from a very 
small number of positively charged trap states according 
to the Fermi-Dirac occupation function. Temperature 
dependent screening effects were not considered in this 
calculation as they give only a change in p by one or two 



orders of magnitude, but the calculated resistivity is by 
about 30 orders of magnitude too small compared with 
the measured values. In the same density interval, the 
measured resistivity decreases just by about one order 
of magnitude, which appears as a nearly flat line on the 
scale of Figure [Io[ 
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FIG. 10. Comparison of p(n) at T; = 3.1 K for experimen- 
tal data (circles) and calculations (solid lines) according to 
the AM-trap model. 

We conclude on the AM-trap model that it cannot ex- 
plain the experimental behavior of the SOI (and similarly 
of the Si-MOS) samples at intermediate temperatures as 
the carrier freeze out on the trap states occurs by far too 
fast if the trap maximum is below Ep. Further improve- 
ments of the AM-model may be possible, but one has to 
note that the addition of trap states at other energies or 
the energetical broadening of the trap states around E t 
destroys the critical behavior around n c and the modified 
model would probably not be able to describe the abrupt 
metal to insulator transition. 

Moreover, case (a) of the trap model by AM seems to 
be able to explain the p(T) dependence only in a very 
narrow density range around n c but can not explain the 
metallic behavior for the wide n range as observed in 
Si-MOS and SOI samples. Thus it seems to be possible 
that the hole trap model in the present form may ex- 
plain some important contributions to p(T) in Si-MOS 
and SOI but is unable to describe the whole temperature 
dependence. As mentioned above, the AM model is espe- 
cially prepared for samples with an Si02 oxide layer as on 
border of the 2D electron gas. The model will thus prob- 
ably not be able to explain the metal-insulator transition 
in other material systems with strong changes in p(T) 
like p-type AlGaAs or Si/SiGe. On the other hand, it is 
not clear anyway whether the apparent metal-insulator 
phase transition has a universal origin or if it is caused by 
several different effects in the different material systems. 
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C. Interband scattering 

The third suggestion (iii) is attributed to interband 
scattering between different conduction or valence bands. 
The correlation between the existence of two conduc- 
tion bands and the metallic behavior was proposed by 
Pudalov for Si-MOS.tfl He assumed that the spin-orbit 
coupling due to the asymmetry of the inversion layer 
potential is strong enough in order to create the split- 
ting. But the calculation of the spin-orbit coupling is 
not straight, forward and have to be performed, with 
great care.L3 A later estimate by Pudalov et al.E£l gives 
a clearly smaller value for the spin-orbit coupling con- 
stant of ol w 6 x 10 -6 K cm, so that the corresponding 
splitting should not be the origin of the metallic behav- 
ior in Si-MOS structures. For p-type AlGaAs Papadakis 
et al. have arguedtH on the spin-orbit coupling due to 
the inversion asymmetry of the electron layer potential, 
but it turned out that the effect is due to an anomalous 
behavior of the magneto-oscillationsH 

Nevertheless, the apparent metallic behavior in p-type 
AlGaAs layers strongly points towards a two-band trans- 
port effectEj There, the heavy hole band is split and 
inelastic interband Coulomb scattering between the two 
bands with different dispersion takes place. This was con- 
firmed by magnetotransport measurements, where the 
positive magnetoresistance is attributed to the classical 
two-band scattering effect. But for n-type Si-based het- 
erostructures, the two-band transport should not be of 
great importance, as the conduction band splitting seems 
to be rather small. 



D. Temperature dependent screening 

The forth suggestion (iv) on temperature dependent 
screening evolves from calculations of the mot 
sus electron density-, dependence for Si-MO£ 
Si/SiGeEJ and III-VEJ structures in the 80'th. In these 
early works, the resulting p(T)-dependence was assumed 
to be relative small. The temperature dependence of the 
resistivity is dominated by the behavior of the screening 
function S(q). For T = 0, this function has a distinct 
kink at q = 2kp with fcp the Fermi-|Wave vector. For 
higher T, this kink is smeared outE3 and changes the 
scattering efficiency especially for q — \k — k'\ « 2k F . 
Thus for large angle scattering, the temperature depen- 
dence of the screening function strongly influences the 
scattering efficiency. 

Gold and Dolgopolov give analytic results for the 
screeningQ that consist of linear-in-T and T 3 / 2 terms 
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Sa(T) w -C(a,n) 



E f 



D(a, n) 



k R T 



-i 3/2 



E f 



(1) 



is dominated by impurities (a = — 1) or by surface- 
roughness (a = 0). We have fitted the p(T) behavior 
with the Eq. [|, assuming a = — f and C(a) — 2 In 2. The 
prefactor C(n) was used as a fitting parameter. As the 
analytic results are only valid for fc^T < Ep/A : the fit 
was limited to that range. Figure |lf| shows that the fit 
gives a very good agreement with the observed T depen- 
dence. 
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FIG. 11. Fit of the temperature dependent resistivity with 
T and T 3/2 screening terms for k B T < Ep/i. 



with C{a,n) = 2C(a)C(n), D(a,n) = 2.45C(a) [C(n)f 
and C(a) a constant depending whether the scattering 



The fitted prefactor C(n) increases from 0.4 at n = 
10 x 10 11 cm~ 2 to 0.55 at 2 x 10 11 cm~ 2 . For lower densi- 
ties, the fit does not give reasonable results as the slope 
of p(T) changes sign at the critical density n c which can 
not be described by the screening approximation alone. 
In Fig. |l2|, the range of the increase in C(n) towards 
lower densities is compared with the-calculated behav- 
ior according to Gold and DolgopolovlD. As can be seen, 
there is good agreement with the curve of the local field 
correction parameter of G(2kp) = 0.75. But if one 
inserts 2k F for q in G(q) = (l/2g ly )g/(g 2 + kj,) 1 / 2 as 
given for Hubbard's approximation in Ref.Q, one obtains 
G(2kp) = 0.224 (with g v the valley degeneracy). Thus 
the above screening description is able to describe the 
temperature dependence but has to use a larger value 
for the local field correction parameter G(2kp) as follows 
from the Hubbard approximation. This discrepancy pos- 
sibly results from the fact that in local field correction 
parameter only exchange effects are taken into account 
and no correlation effects. But the latter might become 
important for large electron-electron interaction parame- 
ters r s . The parameter r s describes the ratio of Coulomb 
to kinetic energy and can be written as r s = g^/dEy/im, 
with as the Bohx. radius. Recently, Das Sarma and 
Hwang calculatedEj that due to temperature dependent 
screening, changes in p(T) by even up to one order in 
magnitude can be explained by numerical calculations. 
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FIG. 12. Density dependence of the prefactor C(n) in the 
screening fit. 

Another puzzling phenomenon in high-mobility Si- 
MOS structures, but also in the SOI samples investigated 
in this work and in other material systems, is the strong 
positive magnetoresistance in parallel magnetic field Bu . 
In Si-MOS structures a more than ten-fold and in SOI a 
nearly four-fold increase in p(Bu) was observed. Dolgo- 
polov and Gold calculatedcJ that an increase of up to a 
factor 4 in p(Bn) can take place as the electron system 
becomes completely spin polarized. This strong increase 
is caused by changes in the density of states at Ep , in the 
Fejapi wave vector Uf and in the screening wave number 
qjEl. We have fitted our data on SOI according to this 
theory, the results are shown in Fig. O. 
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FIG. 13. Fit of the SOI magnetoresistivity p(B\\) with the 
screening behavior according to Dolgopolov and Gold. Den- 
sities are the same as in Fig. |^. The inset shows the density 
dependence of the ratio a = g*m* / gm according to the per- 
formed fit (x-symbols) and to the literature (solid line). 

As can be seen, at high electron densities (i.e. low re- 
sistivity) the fit in Fig. |l^ is in good agreement with the 
data, whereas at lower densities (i.e. high resistivity and 



strong changes in p(B\\)), the screening behavior leads to 
curves which are rather linear and do not show the cur- 
vature of the measured values at low Bu fields. For low 
density, at high Bu fields, the model calculation shows 
an abrupt saturation of the resistivity which corresponds 
to the complete spin polarization of the electron system. 
Both discrepancies between calculation ancLmeasurement 
may be caused by the fact that the modelED was formu- 
lated for T — OK. In that case, the Fermi-Dirac distri- 
bution is a step function and leads to the kink at the 
-B|| saturation field. The approximately linear behavior 
of p(M\\) for small B\\ is probably caused by the distinct 
kinkE3 in the screening function at zero T. For B = 
only one kink exists at q = 2kp which splits into two 
kinks for B > 0.E3 At T > 0, the kinks in the screening 
function are rounded and thus a small splitting of the 
two screening functions (by B|| ) for the spin up and spin 
down electron systems will not be so effective. 

In order to fit the screening model to the data, an ef- 
fective g* factor was used as the only fitting parameter 
instead of a constant value for g. The resulting ratio 
a = g*m*/gm is shown by cross symbols in the inset of 
Fig. h3 in comparison to the ratio a from the analysis ul 
Shubnikov-de Haas measurements on Si-MOS samplesES 
(solid line). Over all, there is a good agreement between 
a from the screening fit and the Si-MOS values. Just 
at the lowest n, the deviation is non-monotonic, but 
this may be caused by the large deviations in the cal- 
culated and measured p(Bu) values as the calculation is 
performed for 0. On the other hand, the above 

screening modelEU is also limited by the Hubbard ap- 
proach, which does not take into account electronic cor- 
relation effects. A further work by Gold indicatescj that 
even larger ratios of p{B\\ > B c )/p(B = 0) than four can 
be explained beyond the random phase approximation 
by taking into account exchange/correlation (i.e. many 
body effects) and multiple scattering effects . In this 
framework possibly the data for high-mobility Si-MOS 
samples can be explained as well. 



E. Ballistic interaction corrections 

The fifth and last suggestion (v) on the origin of 
the metallic state in 2D concerns coherent interaction 
corrections (CIC) to the conductivity in the ballistic 
regime. The theory for it was recently introduced by 
Zala, Narozhny, and Aleiner (ZNA)B. Previously, Alt- 
shuler and Aronov have calculated quantum corrections 
to conductivity (which include especially the interaction 
corrections) in the diffusive regime (i.e. at low temper- 
atures, where ksT <g; h/r) and found a typical m(T) 
behaviorEJ. For the intermediate T range (i.e. the bal- 
listic regime, whece fc^T > h/r), ZNA found a linear 
in T contribution.^! In addition, they calculated the co- 
herent interaction corrections for the transition regime, 
where the behavior changes from diffusive to ballistic. In 







the ballistic regime, the dominant processes can be un- 
derstood as coherent backscattering of electrons at the 
Friedel oscillations o£ the other electrons around some 
perturbation centers.u In contrast, the coherent scatter- 
ing in the diffusive regime occurs at two or more centers 
into a wider range of angles. In the approach of ZNA, 
where the interaction corrections are calculated in terms 
of the Fermi liquid parameter Fq , correlation (Hartree) 
and exchange (Fock) terms are included. As the Fock 
terms have opposite sign to the Hartree terms and the 
relative strength of the two contributions varies with tem- 
perature and electron density, the resulting interaction 
corrections may be either positive or negative. The rela- 
tion of the ballistic interaction correction to the former 
screening approach will be discussed later. 

We have fitted the temperature dependence of the SOI- 
resistivity by the theory of ZNA, but adopting the situa- 
tion for two-fold valley degeneracy in n-type Si inversion 
layers. As was shown by Punnoose and Finkelsteincd, the 
valley degeneracy leads to 15 spin triplet channels instead 
of 3 for the one- valley situation. In both situations there 
is just on charge channel participating. This modifica- 
tion of the ZNA theory was also taken into account in 
the evaluation of the resistivity in Si-MOS by Vitkalov 
et al.ci and Pudalov et alH3. 

In order to fit the experimental data p{T), the func- 
tions /(TV) and t(TT;Fo) for the charge and triplet 
channels, respectively, were calculated according to Ref.El 
and tabulated. The Fermi liquid parameter Fq and 
the conductivity <j(T — 0) are used as fitting parame- 
ters in a true multi-parameter least-square-fit procedure. 
As in the screening case (iv), the fit is performed for 
k B T < Ep/A. It can be seen from Fig. [ll] that the fit- 
ted curve (solid line) describes the experimental values 
(circle symbols) quite well. The inset of Fig. |l4| shows 
the -F^-values as obtained from the fitting procedure. In 
the density range of 2 to 10.4 x 10 11 cm~ 2 the Fq values 
lie between —0.23 and —0.31. For the lowest density of 
1.67 x 10 11 cm -2 , we find a significantly higher value of 
Fq = —0.18. As can be seen from Fig. |l4|, the fit was 
performed only in a short temperature interval and the 
fit is not so significant there. On the other hand, one sees 
also that slop of the T-dependence is indeed drastically 
changing and one eventually comes into a different regime 
of the sample behavior. For even lower electron densities 
the fit did not give any reasonable results as the sample 
finally becomes insulting below n « 1.5 x 10 11 cm -2 . 

In our fitting procedure- we use c(0) as a fitting pa- 
rameter. In the work ZNAB and of Vitkalov et alMd, it is 
suggested to perform a linear extrapolation of the resis- 
tivity p(T) towards zero T. But this has the disadvantage 
that the other contributions (e.g. logarithmic ones) may 
lead to an offset and the value obtained be the extrapo- 
lation towards zero is not the optimal one. As we have 
included the diffusive (i.e. logarithmic) terms in our fit, 
we find a discrepancy between the fitted value cr(0) and 
the one from the linear extrapolation. With our proce- 
dure, we get a better consistency for the fit. 
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FIG. 14. Fit of the temperature dependent resistivity by 
the electron-electron interaction corrections in the ballistic 
regime according to Zala et al. 

If we perform our fit on p(T) under the assumption of 
no valley degeneracy, then the resulting Fq values have 
to be further negative and are typically in the range of 
—0.45 to —0.5. This comes from the fact that the single 
charge channel is always localizing whereas the triplet 
channels contribute the delocalizing parts to the conduc- 
tivity and the strength of the triplet channels increases 
towards negative values of Fq . If there is just one val- 
ley present, the Fq parameter has to be larger in order 
to lead to the experimentally observed delocalizing be- 
havior towards low T whereas in the case of two valleys 
with 15 triplet channels, the Fq value can be lower and 
still leads to the delocalizing behavior. For one valley the 
border between localizing/delocalizing behavior is given 
by Fq w —0.4, for two valleys it is about —0.2. From the 
quality of the fit, it has to mentioned that the screen- 
ing fit of (ii) is nearly identical, the differences are very 
small. If one calculates the prefactor of the linear term 
in the ZNA theory and compares it with C(n) from the 
screening fit, one gets very similar values. This means 
that that in our case the Fock contribution seems to be 
not very important. 

The interaction corrections at intermediate temper- 
atures were also .calculated for an in-plane magnetic 
field Bu by ZNAE Due to the bare Zeeman splitting 
E z = gpsB, the triplet states with different total spin 
component L z are split and the interferenc9-,effects (i.e. 
the triplet channels) are partly suppressed.Ei Note that 
one must not take the renormalized Zeeman splitting 
E* = g* pbB, as homogeneous collective modes are not 
renormalized by electron-electron interaction^ We as- 
sumed again a two-fold valley degeneracy. In addition 
to the effective Zeeman splitting, a possible valley split- 
ting also contributes to the conductivity corrections- The 
expression for the different corrections is given bytiil 

5a ee = Sa c + l5S(j T + 25a(E z ) + 2Sa(A v ) 
+ Sa{E z + A v ) + Scr(E z -A v ), 
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where the terms correspond in the order of appearance 
to the charge and triplet channels, to Zeeman and valley 
splitting and to combinations of the latter two. 

The behavior of the triplet states in parallel field 
and/or valley splitting is described by the functions 
Ki,(x,Fq) and Kd(x/a, Fq) for the ballistic and diffu- 
sive case, respectively^ Here x = E X /2T and E x is any 
of the energies given as an argument in Scr of Eq. ||. We 
have again calculated and tabulated these functions in 
order to perform a true multi-parameter least-square fit 
procedure. First we describe a fit for Ay = 0; the case 
of Ay > is discussed later. For the case of Ay = 0, the 
conductivity corrections reduce to 



5a ee = 5a c + 15&7X + 48a (E z ). 



(3) 



Each Sa(E z ) is able, to freeze two triplet channels 2Sar 
(see Eq. 11 in ZNAQ), so that in the high magnetic field 
limit E z 3> IcbT only 7 triplet channels remain. 



Figure 15 presents the result of the fit with Fq and p(0) 
the fitting parameters and Ay = 0. As can be seen from 
the figure there is a very good agreement for high densi- 
ties over the entire B\\ -field rang and also for the lower 
densities for the low field part. At high parallel magnetic 
field and low density there is a discrepancy between cal- 
culations and experiment. This can be understood, as a 
limit of applicability of Eq. || is that E z < Ep. For the 
lower electron densities in our experiment, we reach the 
total spin polarization regime and the latter relation is 
not fulfilled any more. At first we fitted the whole Bn- 
range, but after recognizing the deviations at high fields, 
we limit the fit range to Bn < 4T as is shown by the 
full lines in Fig. [l5| The extrapolated p(B\\) behavior is 
indicated by the dashed part of the lines. 

The inset of Fig. [l5| depicts the obtained values for the 
fitting parameter Fq , which are in the range of —0.21 
to —0.28 for the higher electron densities. For the lower 
densities, we find a sudden upturn of Fq(u). It occurs 
for the same densities, for which there is a large and sys- 
tematic deviation at high magnetic fields. Possibly, we 
reach a limit of the applicability of the theory. 

We have also performed a least square fit with the 
additional fitting parameter Ay for the valley splitting. 
We get values between ss and 2K for Ay. These val- 
ues do not depend systematically on n, but are rather 
scattered around in the before mentioned range - the 
significance of the parameter Ay is very low. If we set 
Ay = 0, as used for the fit above, the deviation of the 
fitted curve from the data points is nearly the same and 
the fitted values of Fq and p(0) change only marginally. 
This means that we do not get a reasonable estimate 
for Ay from the fit of our magnetoresistivity data for 
in-plane magnetic field. 
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FIG. 15. Fit of the SOI magnetoresistivity p(-B||) with 
electron-electron interaction corrections in the ballistic 
regime according to Zala et al. Densities are the same as 
in Fig. []. The inset shows the Fermi liquid parameter Fq 
(circle symbols) in comparison to the values calculated from 
the effective g-factor in Si-MOS (solid line). 



Figure 16 shows a comparison of the Fq values as ob- 



tained from the fits of the T- and the Bii -dependence. 
In the electron density range of 2 to 6 x 10 11 cm~ 2 , the 
values (marked by circles and squares) are close together. 
This means that two independent measurements and fits 
lead to very similar Fermi liquid parameters Fq. This 
observation supports the validity of the ballistic inter- 
action theory. For n < 2 x 10 11 cm -2 , the Fq val- 
ues suddenly increase strongly. It seems that in this 
regime the description of the experiment by the inter- 
action corrections in the ballistic regime breaks down. 
This is in good agreement with the observation in Si- 
MOS, that for n < 2 x 10 11 cm -2 , the strong decrease 
in p(T) towards lower T is in the diffusive regime, so 
that the ballistic regime is active only for higher elec- 
tron densities O For comparison, Fq was also calculated 
and plotted for .high-mobility Si-MOS structures (solid 
line) as obtainedEa from the effective <?*-factor by the re- 
lation g* = 2/{F£ + 1). The Fq values obtained by the 
two different methods show the same trend in that they 
increase towards negative values for decreasing electron 
density, although the value appears to be more negative 
for Si-MOS. This discrepancy is not clear at the moment. 

If ones compares the fit of our SOI magnetoresistance 
data by the theory on interaction corrections at interme- 
diate temperatures (ZNA)0 with the theory on temper- 
ature dependent screening (Dolgopolov and Gold)E2l, it 
can be seen that the former theory gives a better descrip- 
tion of the experimental data at low magnetic fields. In 
the first case, the nearly parabolic dependence of p(B^ ) 
around B — is described well, whereas in the second 
case the fit gives a nearly linear By-dependence. But as 
discussed earlier in this work, this is not astonishing, as 
the theory on screening in parallel fields is worked out so 
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far for T = 0, where the screening function has the well 
defined kink at q = 2k p . 
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FIG. 16. Comparison of the Fermi liquid parameter Fq 
as obtained from the fit of the temperature dependence (cir- 
cles) and the in-plane magnetic field dependence (squares) 
of the SOI resistivity. 



F. On a possibly new ground state 

The first suggestion (i) on a new quantum mechanical 
ground state in the metallic regime should be discussed 
once more in the light of screening and coherent electron 
transport in the ballistic regime. In a recent workO, we 
concluded that the metallic state in Si-MOS samples is 
not caused by quantum coherent effects. We discussed 
screening as a prominent candidate for the metallic ef- 
fect. At this time screening was not considered as a quan- 
tum coherent effect, primarily. In the new sight of ZNA, 
screening can be described by coherent back scattering 
of electrons at the perturbation centers and it's Friedel 
oscillations. It seems like suddenly the effect has changed 
it's habits. 

If we now try to describe the strong resistivity decrease 
in a combined picture of screening and coherent ballistic 
contributions (let's denote it coherent screening regime), 
one should say that the correlation (Fock) contribution 
may give an important contributi on a gd was missed in 
the former screening descriptions!^!^!! Only in some 
later work, the correlation contributions were included 
in the local field contributiorO'EJ. On the other hand, 
the effect still reflects the behavior of the screening func- 
tion with it's sharp kink at T — which is rounded for 
T > 0. In any case, the coherent screening effects leads to 
a linear T behavior and no new or extraordinary effects 
will emerge from itself. 

The important realization is that the coherent conduc- 
tivity corrections in the diffusive regime can be under- 
stood as a natural continuation of the coherent screen- 
ing. It becomes important as soon as the coherent 
scattering extends over more than one scattering cen- 
ter at a time (which is known since many years). The 
diffusive regime describes the weak localization effect 



and the coherent disptder-induced electron-electron in- 
teraction corrections Hi Only if the latter contributions 
lead to strong and dominant conductivity corrections, 
one should speak of a new quantum mechanical ground 
state. Such_.a scenario was described by Punnoose and 
FinkelsteinEj in the case of two degenerate valleys where 
the delocalizing effect of the 15 triplet channels can 
be significantly enhanced due to renormalization effects. 
But this is not the case for higher electron densities of 
n > 2 x 10 11 cm -2 as was shown in the recent work an. 
the borders for quantum effects in the diffusive regime.Ell 
Only for lower electron densities, such a renormalization 
effect could significantly increase the delocalizing contri- 
butions and lead to a new ground state. 



IV. SUMMARY 



In conclusion, we have shown that silicon-on-insulator 
structures show a very similar 2D-metallic state as Si- 
MOS structures with comparable peak mobility. We dis- 
cussed and compared several models on the physical ori- 
gin of the metallic state. The theory on coherent inter- 
action corrections at intermediate temperatures is able 
to describe the observed temperature and in-plane mag- 
netic field behavior over a wide density range. Under the 
assumption of a two-fold valley degeneracy, the Fermi- 
liquid parameter Ffi is between —0.21 and —0.31. The 
derivated Fermi-liquid parameters from the temperature 
and the in-plane magnetic field dependence are in very 
good agreement to each other. This observation points to 
the validity of the ballistic interaction theory. We do not 
get a significant influence from the valley splitting pa- 
rameter Ay. The earlier screening approach gives also a 
similarly good description for intermediate temperatures. 
This indicates that the correlation (or Fock) terms may 
not be very significant in the investigated regime. The 
charged trap state model, as given in the literature, is 
not able to describe the resistivity versus density behav- 
ior further away from the critical density of the metal-to- 
insulator transition. Also the two-band transport model 
seems to be not applicable for electrons in Si-layers. It 
is not clear yet whether there exists a true metallic state 
for low temperatures at low electron densities or not. 
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APPENDIX A: TRAP-MODEL CALCULATIONS 

Analytic results for the metal-insulator transition in 
gated semiconductors according to trap states-were cal- 
culated recently by Altshuler and Maslov (AM)Ej. In this 
work, AM assume that a hole trap exists in the oxide 
layer, which can easily be charged and discharged. The 
trap state is homogeneously distributed in space but has 
a distinct energetical position. Some of the equations 
used by AM were not derived in that work according to 
the length limits of the letter publication and in addition 
several approximations were used there. We find it there- 
fore useful to discuss the calculations and our numerical 
results in more detail here. 

Figure [T^ shows a simple scheme of the hole trap states. 
If the chemical potential fi (we use here now the same 
notation as by AME3) is below the trap state energy e t 
(Fig. |l~7|a), a hole is trapped, it is positively charged and 
acts as an efficient scattering center. For the case that 
/j, is above £ t , the center is neutral and is inefficient for 
scattering. 
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FIG. 17. Schematic representation of a hole trap center in 
relation to the 2D electron gas. For fi < et the trap state is 
positively charged whereas for /i > et the trap are neutral 
and do not act as scattering centers. 

For the scattering centers, a dipole approximation is 
used by AM. It is assumed that a charge in the ox- 
ide layer will introduce an opposite mirror charge in 
the 2D layer due to screening on in-plane electric fields. 
For the energetical position of the trap state, a term 
e 2 /(87reo£ox^) has to be taken into account, with eo = 
8.854 x 10 12 C 2 /Nm 2 the vacuum permeability, e ox the 
relative dielectric constant of the oxide and z the distance 
of the charge from the 2D electron layer; the calculations 
are performed in SI units here. The energetical position 
of the charged trap states also depends on the applied 
gate voltage via eVgz/d, with d the thickness of the ox- 
ide layer. Figure [18| depicts the dependence of the total 
energy of the charged trap states in the dipole approxi- 
mation as a function of the distance z for different gate 
voltages V g . Only in the spacial region where et(z) > /i, 
the centers are positively charged. 

The expression for the resistivity p follows in the 
Drude-Boltzmann approximation by calculating the ef- 
fective transport scattering time 



J der(e)edf/de 
J deedf/de 



for the conductivity a — ne 2 (r) / m* . The weight- 
ing of (t) with the kinetic electron energy e in Eq. 



Al fundamentally follows from the Drude-Boltzmann 



approximation!.-!, as the Fermi velocity vf and the shift 
of the Fermi surface in k-space are both proportional 
to yfe, which enter into the expression for the current 
j x — — e J dkn(k) f (s)v x . The integral in the denomina- 
tor of Eq. ^l] is in 2D just equal to the Fermi energy Ep 
- also for elevated temperatures. 

trap energies in oxide for fixed Etc-EO 




FIG. 18. Full lines show trap energy vs. distance z from 
the 2D electron gas layer plotted for different gate voltages 
V g . The dashed lines indicate the potential energy of the 
trap without the mirror charge term; they merge on the right 
hand side at the trap energy eto- The thick horizontal line 
indicates the position of the Fermi energy Ef\ it actually 
consists of several thin lines for different values of _Ef(V' 9 ), 
but as the Fermi energies are so small in comparison to the 
considered trap energies, they are not resolved. The stars de- 
pict the maxima of the energy trap curves, the triangles show 
the mean z-value of the positive region, i.e. where et > £o 
with Eo the lowest 2D energy level in the Si inversion layer. 

The transport scattering time r(e) has to be calculated 
from the scattering rate 



f d 

!/ T ( e ) = / dzN^(z)v(e)<j scm (e,z) 
Jo 



(A2) 



with iV t 1j(z) the three dimensional density of charged 
traps, v(e) the electron velocity and <7 scm (e,2:) the 
scattering cross section for momentum relaxation. Ac- 
cording to AM, the scattering cross section for a dipol 
is given by cr scm (e,z) = 2.74(e 2 z 2 /87re e*e) 1 / 3 . The ef- 
fective relative permitivity e* = (esi + e ox )/2 ~ 7.9 is 
the mean value of the Si and the SiC>2 layer. The den- 
sity N^(z) = N t3 P + (z) in the AM notation with N t3 
a constant density of traps per volume and P+{z) — 
[0.5 exp((fi — e e (z) — et) / ksT) + , the thermodynamic 
probability for a trap state to be chargedEj. 

By inserting the above expressions into Eq. |A2| one gets 



(Al) 



l/r(e) = c'N+zT^ 6 



(A3) 
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with the factor d = 2.74(e 2 /87re e*) 1 / 3 ^/2/rn*, an ef- 
fective number of positive trap states per area ATg = 
N t 3 (Az) cS , the effective width of positive charge layer 
(Az) eff = J dzP + (z)(z/z m ) 2 / 3 and the position z m = 
y^ed/S^£oeoxVg of the energetical maximum of the trap 
energyEI By inserting Eq. A3 into Eq. Al , one gets 



(r) oc 1/fi / dee 5/6 df/de 



(A4) 



For the Drude resistivity p — m*/ne (t), an effective 
energy e can be defined so that a relation equivalent 
to Eq. A3 can be written for the effective values, i.e. 
1/ (r) = c'N^Zm 3 ^ 1 ^ 6 - A simple calculation gives 



e = p 



dee 5/6 df/de 



(A5) 



The first derivative of the Fermi-Dirac function / can be 
expressed as 



(A6) 



df/de = - /(l - /) 

= - (4fc s Tcosh 2 ((£ - f i)/2k B T))- 1 

and one obtains the same expression as Eq. 8 in Ref.E3 
With these relations, the resistivity can exactly be writ- 
ten in terms of the effective energy s as 



p = — ~N^v{e)cr scm {e, z n 



(A7) 



We have calculate the resistivity p(n, T) by numeri- 
cal integration for the effective thickness Az of the layer 
with charged-trap states and the effective energy e. In the 
work of AME3, a parabolic approximation was used for 
the dependence of the trap energy e t on z and an analytic 
expansion was given for the behavior of p{n,T). With 
this, AM were able to show that in principle a critical 
behavior in the density dependence of p(T) can appear. 
As in our work, we want to test whether the metallic 
behavior can be explained for a wide density range or 
not, we use the numerical calculation method in order to 
be able to perform a precise evaluation of the model be- 
havior for large deviations of n from the critical density 
n c . 

AM proposed in their model that the dependence of 
p(n) in the vicinity of n c comes from the strong change 
of the number of positively charged centers, as the Fermi 
energy is just above or below the energy maximum e m 
of the trap states, where they have a singularity in the 
density of states. The sensitivity on temperature comes 
from the broadening of the thermodynamic probability 
distribution P + with increasing temperature. In order to 
test this behavior for the strong temperature dependence 
of the resistivity, i.e. the most characteristic feature of 
the metallic state, we calculated the temperature depen- 
dence within the model at an intermediate temperature 
of T = 3K. For densities n < 5 x 10 11 cm -2 , p(T) does 



not saturate yet at its low-T value. The saturation be- 
havior is not described by the AM-model. The low-T 
value has to be cause by other mechanisms like residual 
impurities which do not change their charge state in the 
same manner. For the evaluated value of T 3K, the 
T-dependence should behave according to the AM-model 
if it is able to describe the metallic state. 

The result of our numerical calculation is shown in 
Fig. 
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It is clearly visible there, that the calculated p- 
dependence changes by many more orders of magnitude 
than the measured one. In order to analyze the origin 
of the strong change in the calculated p{n) behavior, we 
show in Fig. [1^ the dependence of the effective number 
of charged trap states per area N^~ s on the electron den- 
sity (i.e. the gate voltage V g ). As long as e m > p a part 
of the trap states is positively charged even for T = 
and the variation in effective positive charge layer width 
Az is small. But for e m < p, the trap states are only 
charged according to the exponentially decaying part of 
the thermodynamic probability function P+(z) and thus 
the strong exponential T-dependence occurs. 

In the same density range, the other v alue s which enter 
into the resistivity calculation (see Eq. A7) change only 
by about a factor of two. It is thus clear that within the 
proposed AM model, the strong shift of the energetical 
position of the trap states with the applied gate voltages 
leads to too strong changes in p(n) further away from the 
critical density n c . 




n (10 cm ) 

FIG. 19. Effective number of charged trap states per area 
vs. electrons density (i.e. varying gate voltage) at tem- 
peratures of T = 1, 3, and 5 K. 
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